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Abstract 

Does quantum dynamics play a role in DNA replication? What type of tests would reveal that? Some statistical checks 
that distinguish classical and quantum dynamics in DNA replication are proposed. 
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Genetic information processing 

Darwinian evolution — survival of the fittest — is an opti- 
misation problem. Its foundation is the observation that 
living organisms have adapted to their environment, and 
have exploited the available material resources and the 
physical laws governing them to the best of their capabil- 
ity. The environmental influences are usually local and 
so are the adaptations. So the meaning of "best of their 
capability", is different in different situations, and it is 
an interesting exercise to figure out the physical princi- 
ples involved in the optimisations. 

The genetic code is a feature of living organisms that 
has remained virtually unchanged from the ancient bac- 
teria to modern human beings. Its purpose is to convey 
the hereditary genetic information from one generation 
to the next. The task is carried out by the long chains of 
DNA molecules residing in the nuclei of the cells, which 
faithfully replicate themselves each time the cell divides. 
In the early studies of molecular biology, it was proposed 
that the genetic code was a frozen accident (Crick 1968). 
More recent investigations have shown that some parts 
of the genetic code have flipped back and forth indepen- 
dently in completely different organisms (Knight et al. 
2001). This leads credence to the belief that the genetic 
code is evolvable and optimisation criteria have been in- 
volved in its design. 

What are then the optimisation criteria involved in 
genetic information processing? The task involved is an 
unsorted assembly operation: pick up the building blocks 
from their random mixture in the environment, check 
whether they have the desired properties or not and then 
join them together in the requisite order. This is a well- 
defined computational problem, and it is reasonable to 
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suppose that optimisation corresponds to carrying out 
the task as quickly as possible and without errors. The 
solution to this problem depends on the type of dynam- 
ics used to implement the algorithm; in particular the 
optimal solution is different for classical and quantum 
dynamics. 

I have proposed that assumption of quantum dynam- 
ics in the assembly operation would naturally explain 
why the genetic information is organised in a base-4 lan- 
guage, and not in the base-2 language of our classical dig- 
ital computers (Patel 2000b). My study of the molecular 
structure of DNA showed that it indeed has the capabil- 
ity to implement the quantum algorithm (Patel 2000a). 
The question that remains is: does DNA replication ac- 
tually use the quantum algorithm, or is the existence of 
both software and hardware features an accident? 

Of course, the structure of DNA came into existence 
billions of years ago, and it could be that what was rele- 
vant when life arose is not relevant now; i.e. the observed 
features are hanging around from a bygone era like the 
human appendix. Unfortunately, we cannot recreate the 
conditions in which life originated, and so cannot test 
whether quantum dynamics played a part in it or not. 
What we can test is whether quantum dynamics plays 
a role in DNA replication today or not. Direct observa- 
tion of the assembly operation as it happens would be 
an ideal test, but our technology has not yet progressed 
to that level. We have to rely on statistical tests to infer 
the dynamical process. 

Here I point out some statistical tests, that can dif- 
ferentiate between classical and quantum dynamics in 
DNA replication. They crucially depend on our ability 
to construct "designer DNA" and use it in the replica- 
tion process. The following section describes tests based 



on DNA replication rates, while the subsequent section 
describes tests based on the DNA structure. 

Replication rate checks 

First let us go through the mathematical description of 
the assembly algorithm. That clearly brings out the dis- 
tinction between classical and quantum dynamics, and 
also explains the criterion for optimisation. 

Classical unsorted assembly process 

The assembly process is a variation of the search process, 
where the objective is not to locate the desired item in 
a database but to assemble it using the available build- 
ing blocks. The building blocks can be of many different 
types distinguishable from each other by certain proper- 
ties. Let a be the number of types of building blocks and 
n be the length of the chain to be assembled. Then the 
total number of items that can be assembled is N = a n . 
N is the measure of information contained in the chain, 
and a > 2 is necessary to be able to convey worthwhile 
information. 

We consider the situation where building blocks are 
available in large numbers for each type. Also they are 
available not in a particular order, but as a random en- 
semble of all possibilities. Let td be the time required to 
pick up a building block from the random ensemble, and 
t a be the time required to add each building block to the 
growing chain. Since the ensemble of the building blocks 
is random, the building block that is picked up may or 
may not be the one desired for adding to the growing 
chain, and therefore td 7^ t a . Certain property checks 
decide whether the building block picked up is the de- 
sired one or not. If it is the desired one, it is added to the 
growing chain. If it is not the desired one, it is discarded 
and another one is picked up from the random ensemble. 
Again its properties are checked, and the trial and error 
process continues till the desired building block is found. 
We also assume that the property check process has no 
memory, i.e. there is no correlation between a building 
block that is discarded and the new one that is picked 
up next. The ensemble of building blocks is thus truly 
random; what is discarded once can be picked up imme- 
diately again. In this situation, the probability that the 
picked building block is the desired one is 1/a. 

The total time required to classically assemble the 
complete chain is (N is chosen to be sufficiently large 
so that round-off errors can be ignored): 

T c (a) = t a n = t d an = t d a\og a N = t d a(\n N/ In a) . (1) 

The optimisation criterion for the assembly process is 
to minimise T c {a) by varying a, while holding td and N 
fixed. Holding N fixed means that the information to 
be conveyed is fixed. Holding td fixed means that the 
total number of all the available building blocks in the 
random ensemble is fixed; the concentration of a par- 
ticular type of building block will be proportional to 



I/a. This minimisation criterion decides the number of 
types of building blocks required for the quickest pro- 
cessing of a fixed amount of information. The best value 
for a obtained by minimising T c (a) is not an integer, 
Omin = e = 2.71828 . . .. The integer solutions close to 
the lowest value of T c (a) are: 

T c (a = 2) = ^tdlnN w 2.8854i d lnAT , 
T c (a = 3) = XtahiNm 2.7307t d lnA , 
T c (a = A) = ^trflnJVw 2.8854t d lnA . (2) 

For larger a, T c (a) increases monotonically. Thus the 
best choice is a = 3, while a = 2 and a — 4 would serve 
equally well as the next best possibility. 

Classical DNA replication rates 

In the particular case of DNA replication, the build- 
ing blocks are the nucleotide bases, and property checks 
are performed by molecular bonds involved in the base- 
pairing, td is the sum of the time required for a base to 
diffuse from the environment into the enzyme cavity and 
the time required for making property checks. The repli- 
cation time per base-pair, t a , can be easily measured in 
experiments. Measurement of td would require observ- 
ing the accept/reject process for a picked base, which is 
not yet directly observable. The analysis of the previous 
section predicts specific DNA replication rates depend- 
ing on the value of a. 

Arbitrary choices of a are not possible with the natu- 
rally occurring molecules, yet situations can be selected 
to correspond to some of the values of a: 

• Genetic DNA of living organisms has a = 4, consisting 
of two complementary base pairs, A-T and C-G. 

• DNA constructed with one of the complementary base 
pairs and a single base on one strand would correspond 
to a = 3. The single base should either pair with itself or 
the environment should contain its complementary base. 
(For example, the strand to be replicated could be made 
of only A,C,G while the environment has only T,C,G.) 

• DNA constructed with only one of the complementary 
base pairs would correspond to a = 2. A — T rich and 
C — G rich DNA already exist in nature, and they can 
be modified to eliminate one of the complementary base 
pairs. 

• DNA constructed with a single base (as in case of 
a = 3) on one strand would correspond to a = 1. Such 
DNA will not carry information, but its replication rate 
can still be measured. 

The classical base-pairing rates as a function of a are 
predicted to be 

R c (a) = n/T c (a) = (l/a)^ 1 . (3) 

The unknown td can be eliminated by forming ratios of 
replication rates. Such ratios can be determined more 
accurately than the absolute rates, since they cancel un- 
known systematic errors in td- For example, it is unlikely 
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that every time the correct nucleotide base diffuses in 
the enzyme cavity, it will bind with the intact strand. 
Wasted encounters with the correct nucleotide base can 
be accounted for by just rescaling td- 

Another possibility for a systematic effect is that the 
times for A— T and C—G base-pairing are different, since 
they involve different free energies. This can be tested 
by comparing the replication rate of a = 2 A — T strands 
with that for a = 2 C—G strands. If there is a systematic 
effect, the replication rates can be renormaliscd to reflect 
the proportion of each base-pair. 

Thus, by constructing DNAs corresponding to differ- 
ent values of a and comparing their replication rates, we 
can experimentally test whether DNA replication follows 
classical dynamics or not. 

Quantum DNA replication rates 

The quantum search algorithm works faster than its clas- 
sical analogue, by using clever interference of a super- 
position of states (instead of individual states) (Grover 
1996). The algorithm can be described as a unitary evo- 
lution of a set of states (one state corresponding to each 
possibility) in the quantum Hilbert space. At the begin- 
ning, any starting state is converted into a uniform su- 
perposition of all the possibilities. Then two steps are ex- 
ecuted iteratively: (i) the amplitude of the desired state 
is flipped in sign by the property checks, and (ii) all the 
amplitudes are reflected about their average value by an 
overrelaxation procedure. These iterations amplify the 
amplitude of the desired state, and diminish the ampli- 
tudes of the remaining states. By stopping the iterations 
at the right moment, the desired state is selected with 
high probability. The number of iterations required, Q, 
is related to the number of items to be distinguished, a, 

by 

(2Q + l)sin- 1 (l/V^) = tt/2 • (4) 

For a given value of a, the solution for Q is not necessarily 
an integer (the exception is Q = 1 for a = 4). For non- 
integral Q, the number of iterations is to be interpreted 
as the integer closest to and larger than Q. 

Let t r be the time required to add one building block 
to the chain using the quantum algorithm. In addition 
to the time required to pick up a building block from the 
environment and performing property checks, it includes 
the relaxation time required for creating the necessary 
superposition state. Thus t r is expected to be larger than 
the time td of the classical algorithm. The total time 
required to assemble the complete chain by the quantum 
algorithm is proportional to the number of iterations, 
and is not a smooth function of a. The lowest values are 
given by: 

T q (a<A) = t r n = t r (lnN/lna) , 
T g (5<a<10) = 2t r n = 2i r (lnA/lna) . (5) 

The minimum value of T q (a) is obtained for a = 4. For 
the DNAs constructed to correspond to different values 



of a, as described above, the quantum base-pairing rate 
is 

R q (a < 4) = n/T q (a) = t; 1 . (6) 

This rate is independent of a distinctly different from 
the classical prediction. Therefore, experiments measur- 
ing DNA replication rates with different values of a can 
distinguish classical and quantum dynamics. 

In reality, the quantum algorithm is unlikely to work 
with 100% efficiency. Disturbances from the environ- 
ment, i.e. decoherence effects, cannot be wished away. 
It should be noted that, with different quantum states 
represented by different physical objects, the quantum 
algorithm needs only superposition and interference, and 
does not need the far more fragile property of quantum 
entanglement (Lloyd 2000). Whenever the quantum al- 
gorithm is disrupted by decoherence effects, the replica- 
tion process can still continue with multiple attempts as 
in the classical case. Even an imperfect quantum algo- 
rithm will beat the best classical alternative, if its success 
probability obeys 

p>T q (a = 4)/T c (o = 3) = 0.264t r /t d . (7) 

This is a substantial tolerance margin for the quantum 
algorithm, provided that t r is not significantly larger 
than td- Biological systems are known to sense small 
differences in population growth rates (e.g. growth of 
type dP/dt = rP). Even an advantage of a fraction of 
a percent is sufficient for one species to overwhelm an- 
other, provided there is long enough time for evolution. 
Nonetheless, it should be noted that an imperfect quan- 
tum algorithm will be harder to verify, because its prop- 
erties will lie somewhere in between the classical and the 
quantum predictions. 

Structural checks 

Quantum superposition of states is a crucial ingredient 
of any quantum algorithm that outperforms its classi- 
cal counterpart. Processing of a linear superposition of 
states is a type of parallel processing which can speed 
up certain algorithms without requiring additional re- 
sources. The superposition of states has to be created 
from the classical state used to initialise the algorithm, 
and thereafter quantum coherence has to be maintained 
during the course of the algorithm. 

The quantum assembly algorithm requires that a uni- 
form superposition of all possibilities be created from 
any classical starting state. This is possible only if the 
quantum dynamics provides transition matrix elements 
amongst the classical states. (The transition matrix ele- 
ments produce a mixture of states in a random classical 
environment, but they produce superposed states in a co- 
herent quantum environment. I have hypothesised that 
polymerase enzymes help carry out these transition pro- 
cesses (Patel 2000a).) The classical states are nucleotide 
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bases with different chemical compositions, so the quan- 
tum transition processes must involve exchange of chem- 
ical groups. Such exchange of chemical groups can be de- 
tected by isotopic tagging; radioactive processes or NMR 
techniques can track the shift of isotopes from one loca- 
tion to another. 

The nucleotide bases consist of a large group of atoms 
common to all of them, and a small group of atoms 
which tell them apart. The superposition state can be 
created by many different types of quantum transition 
processes. But any such quantum transition process will 
need to exchange the small group of atoms between dif- 
ferent nucleotide bases or between a nucleotide base and 
the polymerase enzyme. Consider a situation where one 
of the nucleotide bases in the environment has some of 
the atoms in both the large and small group of atoms re- 
placed by different isotopes. If there are no exchanges of 
chemical groups, then the two isotopically tagged groups 
of atoms will be always found together. On the other 
hand, if the exchanges of chemical groups do take place, 
then the two isotopically tagged groups will get sepa- 
rated, at least in some of the molecules. In such a case, 
after the replication is complete, the isotopically tagged 
small group of atoms will be found either in another 
nucleotide base (most likely incorporated in the newly 
constructed DNA strand) or in the polymerase enzyme. 
This can be experimentally detected. 

Another point worth noting is that if the polymerase 
enzyme plays an active role in the exchange of chemical 
groups, then it is likely to have all the necessary chemical 
groups in close vicinity of its active sites, so that the ex- 
change may take place smoothly. Detailed investigations 
of the polymerase enzyme structure can confirm this. 

Summary 

To take quantum processing of genetic information be- 
yond an explanation, its underlying hypotheses must 
be checked against falsifiable predictions. It is not yet 
technologically feasible to observe the dynamics of DNA 
replication, so I have pointed out some statistical tests 
that predict different results for classical and quantum 
dynamics. The predictions are based on the dynamical 
steps involved in the algorithm, and so they indirectly 
verify the dynamical process. 

The tests are of two types: some measure the repli- 
cation rates and some check the structure of the par- 
ticipating molecules. The replication rate tests require 
the number of building blocks in the DNA strand to be 
variable. The protein coding sections of naturally oc- 
curring DNA have all the four nucleotide bases. On the 
other hand, stretches of DNA strands made up of fewer 
bases are known to exist in the non-coding sections of 
DNA, and it should be possible to extract strands having 
one, two or three building blocks from them. The struc- 
tural checks require introduction of isotopically tagged 
nucleotide bases in the environment, and then tracking 



the various isotopes to their final destinations. Such tests 
have been carried out in many chemical processes. 

These tests use only naturally available molecules, and 
so they check the optimality of the DNA replication pro- 
cess only in a limited sense. In particular, the building 
blocks of DNA are varied, but the enzymes are left un- 
changed. The presence of the enzymes is a must for the 
DNA replication process; random molecular collisions do 
not replicate DNA. It is quite likely that properties of the 
enzymes and the building blocks are correlated, which 
implies that the optimal replication algorithm should be 
searched for by simultaneously changing the enzymes as 
well as the building blocks. I have no clue regarding how 
to find the optimal enzyme. In this sense, the suggested 
tests do not really confirm the optimality of the DNA 
replication process, although the quantum assembly al- 
gorithm was proposed on the basis of optimality. 

Finally, I want to point out that I have used the lan- 
guage of DNA replication for describing the tests, but the 
tests can equally well be carried out for the mRNA tran- 
scription process. The base-pairing algorithmic steps are 
the same in both cases, even though there are chemical 
differences in the building blocks and the polymerase 
enzymes used. In fact, the tests should be easier to im- 
plement for the mRNA transcription process, since it 
involves only one DNA strand (in contrast to both the 
strands being simultaneously involved in the DNA repli- 
cation process). 
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